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Detailed Monte-Carlo simulations were used to evaluate the performance of a liquid neon scintillation detector 
for dark matter and low-energy solar neutrino interactions. A maximum-likelihood event vertex fitter including 
PMT time information was developed, which significantly improves position resolution over spatial-only algo- 
rithms, and substantially decreases the required detector size and achievable analysis energy threshold. The 
ultimate sensitivity to WIMP dark matter and the pp flux uncertainty are evaluated as a function of detector size. 
The dependence on the neon scintillation and PMT properties are evaluated. A 300 cm radius detector would 
allow a ~13 keV threshold, a pp flux uncertainty of ~1%, and limits on the spin-independent WIMP-nucleon 
cross-section of ~ 10~ 46 cm 2 for a 100 GeV WIMP, using commercially available PMTs. Detector response 
calibration and background requirements for a precision pp measurement are defined. Internal radioactivity re- 
quirements for uranium, thorium, and krypton are specified, and it is shown that the PMT data could be used for 
an m-situ calibration of the troublesome 85 Kr. A set of measurements of neon scintillation properties and PMT 
characteristics are outlined which will be needed in order to evaluate feasibility and fully optimize the design of 
a neon-based detector. 
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This paper presents details on the design of 
a dual-purpose liquid neon based detector for 
dark matter and low-energy solar neutrino inter- 
actions. The origin of dark matter in our universe 
and study of low-energy solar neutrinos are two of 
the foremost topics in particle astrophysics. Cur- 
rent observational evidence strongly suggests that 
a large fraction of the total matter in our universe 
is non-luminous, non-baryonic matter outside of 
the current standard model of particle physics. A 
currently favored hypothesis is that the dark mat- 
ter consists of WIMPs (Weakly Interacting Mas- 
sive Particles) . A natural candidate for this parti- 
cle is the lightest superpartner which arises from 
supersymmetry, the neutralino. Current direct 
dark matter searches attempt to detect WIMP 
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interactions through elastic scattering on target 
nuclei, and this is the basis for WIMP detection 
in this study. 

The current generation of solar neutrino exper- 
iments have shown that mixing occurs between 
the neutrino flavors. One of the goals for next- 
generation experiments is to measure with high- 
precision the low-energy pp neutrinos. These neu- 
trinos account for greater than 90% of all solar 
neutrinos, and are well-constrained by the mea- 
sured solar luminosity. A precision measurement 
of the flux of pp neutrinos would lead to improved 
knowledge of neutrino mixing in the solar sector 
(#12), tests of the neutrino oscillation prediction 
in the low-energy region where vacuum mixing 
dominates, and precision tests of solar evolution 
models. 

WIMP particles (x) and solar neutrinos (v) can 
be observed in liquid neon through their elastic 
scattering from nuclei and electrons, respectively: 

X + Nc — ► x + Ne' (1) 
v + e~ — > v'+e~'. (2) 

Monte-Carlo simulations were performed to eval- 
uate the performance of a neon scintillation de- 
tector sensitive to the recoiling nuclei (Ne ) and 
electrons (e~ ) from equations 2] and [3 The con- 
cept of a liquid neon scintillation detector for so- 
lar neutrinos was first proposed by McKinsey and 
Doyle 1 . The experimental design of a next- 
generation WIMP dark matter experiment should 
allow a sensitivity increase of several orders of 
magnitude over present dark matter experiments 
to probe interesting regions of supersymmetry pa- 
rameter space. Ultra-low backgrounds and a large 
fiducial target are required to achieve this goal. A 
large spherical volume of neon, along with precise 
event position reconstruction allows for a large 
fiducial mass and a large mass to surface area ra- 
tio, with the advantage that external source back- 
grounds are physically distant from the fiducial 
target and can be dealt with through position re- 
construction. 

Simulations were performed assuming nomi- 
nal characteristics of scintillation light in liquid 
neon, and detector performance was estimated 
for commercially available photomultiplier tubes 



(PMTs) with ultra-low background glass. The 
details of the simulation inputs and algorithms 
are presented in section 2. The detector re- 
sponse to scintillation events is characterized in 
section 3, where reconstruction algorithms and 
response functions useful for evaluating predicted 
observables for dark matter and solar neutrino 
interactions are presented. The sensitivity to so- 
lar neutrinos is discussed in section 4. Section 5 
describes the reduction of external and internal 
backgrounds through reconstruction and analysis 
cuts, and the ultimate background requirements. 
Optimizing the detector for dark matter and solar 
neutrino sensitivity is discussed in section 6. Cali- 
bration requirements for the full-scale experiment 
are outlined in section 7. The requirements from 
photon absorption and sensitivity to neon scin- 
tillation and detector parameters are presented 
in section 8, along with a list of which measure- 
ments are needed to fully optimize the detector 
design. Conclusions from this work, and poten- 
tial advantages of this approach are presented in 
section 9. 

2. Simulation Details 

A detailed simulation of a liquid neon detector 
was coded using the GEANT4 package 0. This 
package is used for tracking of electrons, 7 and 
x-ray photons, a-particles, recoiling nuclei, and 
neutrons in the detector, and subsequent simu- 
lation and tracking of scintillation photons. A 
3D model of the PMT is used, and tracking of 
photons across surfaces is accomplished with the 
UNIFIED model of the DETECT code 0, as in- 
corporated in GEANT4. All neutrino or dark 
matter interactions arc detected via scintillation 
light produced in the liquid neon. Tabled shows 
the properties of liquid neon used in the simula- 
tion, along with the references for these values. 
Scintillation light is produced via excitation of 
molecular states in the neon. The relative inten- 
sity of these states, as well as the overall yield (or 
quenching factor) depends on the ionizing parti- 
cle. Two states are produced: the singlet ( 1 E) 
state with an assumed 2.2 ns lifetime 0] and the 
triplet ( 3 E) state with a 2.9 ps lifetime 0- 

Quenching factors and the ratio of intensity of 
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the singlet to triplet states ( 1 // 3 /), which de- 
termine the overall time structure of the scin- 
tillation light, for electrons, alpha particles, and 
nuclear recoils were inferred from a combination 
of measurements with liquid argon, xenon, and 
neon |til7| and are shown in Table [21 The dif- 
ference of these intensity ratios between electron 
and nuclear recoils provides discrimination be- 
tween these event types and ultimately allows the 
separation of potential WIMP events from solar 
neutrino events and backgrounds. Many of the 
entries in Tables 1 and 2 are estimates or prelim- 
inary measurements, and need to be quantified 
more precisely. The sensitivity to these parame- 
ters is discussed in section ing To detect the scin- 
tillation photons with a 75% photocathode cov- 
erage on a sphere with a nominal radius of R psU p 
= 300 cm (Rpsup denotes the radius of the photo- 
multipler support structure), 1832 12 cm hemi- 
spherical PMTs are arranged in a 3-frequency 
icosahedral pattern . The nominal PMT prop- 
erties assumed in the simulation are shown in Ta- 
ble|21 The commercially available Electron Tubes 
D737KB 9 PMTs assumed for the simulation 
are made with ultra-low background glass, with 
quoted upper limits of uranium and thorium of 
30 parts per billion (ppb), and of potassium of 
60 parts per million (ppm). These PMTs can be 
manufactured with low-resistivity photocathodes 
for operation at low temperatures |10| . The ran- 
dom count rate in the PMTs is assumed to be 
500 Hz, extrapolated from measurements made 
by the manufacturer . The spectrum of light 
emitted from scintillation in neon peaks at 77 nm, 
and is shown in Fig^a), adapted from |12|. The 
Rayleigh scattering length at this wavelength in 
liquid neon has been calculated to be 60 cm |13| . 
Fig. IHb) shows the Rayleigh scattering length 
versus photon wavelength in liquid neon used in 
the simulation 

The nominal detector geometry used in the 
simulation is shown in Fig. [21 A sphere of PMTs 
is immersed in a volume of liquid neon. Each 
PMT views the inner region of neon, and consists 
of a hemispherical photocathode and glass win- 
dow which has been covered with a wavelength 
shifting film. Extreme ultraviolet (EUV) photons 
which strike the front of the PMT are shifted up 



in wavelength and are emitted into 47r. These 
photons can then be transported through the win- 
dow and generate a photoelectron in the cathode, 
according to the PMT quantum efficiency (0.20) 
used in the calculation. The wavelength shifter 
(WS) is assumed to be transparent to photons 
which have been previously shifted into the 250- 
500 nm range. 



Table 2 



Quenching factors and singlet /triplet intensity ra- 
tios ( 1 // 3 /) for electrons, a's, and nuclear recoils 
(n.r.) used in the simulation. 



Property 


e 


a n.r. 


quenching (prompt) 


1.0 


0.7973 0.25 


i// 3 / 


2.0 


8.7 20 



Table 3 

Nominal PMT properties assumed for simulation 
for Rpsup = 300 cm. 



Property 


Value 


number of PMTs 


1832 


photo cathode radius 


12 cm 


effective coverage 


0.75 


Q.E. (blue photons on envelope) 


0.20 


threshold 


1/4 p.e. 


photo cathode reflectance 


0.02 


dark noise rate 


500 Hz 


Electron Tubes D737KB 


mass (glass) 


1 kg 


thickness (glass) 


3 mm 


uranium 


30 ppb 


thorium 


30 ppb 


potassium 


60 ppm 
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Table 1 

Properties of liquid neon used in the simulation. 



Quantity 


Value 


Reference 


LNe scintillation yield (prompt) 


10791.7 photons/MeV 


|7I14I15I16| 


LNe scintillation spectrum 


(see Fig. 


m 


Rayleigh scattering length @ 80 nm 


60 cm 




refractive index in LNe @ 80 nm 


1.233 


na 


absorption length in LNe @ 80 nm 


CO 




singlet (prompt) time constant (ts) 


2.2 ns 


4 


triplet (slow) time constant (tt) 


2.9 (is 


EH 


prompt WS efficiency for EUV photons (47r) 


1.0 


El 


total WS efficiency for EUV photons (47r) 


1.35 


[T7j 


WS time constant 1 


14.816 ns 




1 Time constant for 100% efficiency within 20 ns, inferred from [7]. 





3. Detector response 



isotropically into 4ir: 



In this section the detector response functions 
in position and energy are evaluated using the 
Monte-Carlo calculations. These response func- 
tions can then be used to estimate the signal 
from WIMP and neutrino interactions, the back- 
grounds present in the signal region, and the ul- 
timate experimental sensitivity. 



3.1. Event position and energy reconstruc- 
tion 

Coakley and McKinsey have shown the poten- 
tial for using completely geometrical position re- 
construction algorithms for scintillation events in 
a liquid neon detector ^Hj- Here we describe 
a more general reconstruction algorithm which 
incorporates PMT hit time as well as hit pat- 
tern information to improve position resolution. 
Although the PMT timing distribution is quite 
broad (many tens of ns) , the distinctive shape of 
this distribution for events near the center versus 
that near the PMT region improves the position 
resolution significantly. 

In the absence of scattering or noise hits, the 
location of a scintillation event can be estimated 
by taking the geometric mean of the PMT hit 
positions, since scintillation photons are emitted 



R 



■event 



N P 
G 



AIT 



NpMT 



PMT 



i=0 



where R even t is the event location and R PMT is 
the vector to the i th PMT to detect a scintillation 
photon. Assuming spherical detector symmetry, 
the only effect of scattering can be a radial offset 
between the true event position and the mean of 
the PMT hit positions, so that we have 



Re. 



where 



= G 



Revent ^ G 
\Revent \ \G\ 



(3) 



(4) 



and a is related to the effective scattering length 
in the detector. Rearranging equation^ we have 



Revent = G(l + a), 



(5) 



absorbing an arbitrary constant and with a in 
general a function of event radius and energy. We 
can thus determine a by Monte-Carlo and use 
equation[S]to estimate the position of scintillation 
events. 

The reconstructed energy and position distri- 
butions (using the spatial algorithm) are shown 
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Figure 1. (a) Wavelength spectrum of neon scin- 
tillation light, (b) Rayleigh scattering length 
in neon. The Rayleigh scattering length versus 
wavelength has been extrapolated from the cal- 
culated value at 80 nm based on the A -4 depen- 
dence. 



in Fig. The reconstruction calibration func- 
tion (a from equation is shown in Fig. 01 The 
position resolution, derived as a function of the 
electron energy, is shown in Fig. [5] where it is 
seen that the position resolution worsens sharply 
as the number of noise hits becomes comparable 
to the number of scintillation hits. 

The energy for each event is estimated by ap- 
plying the correction 

T cS = a + (tix NHITs (6) 

where NHITs is the total number of PMT hits, 



Figure 2. Detector geometry used in the simu- 
lation. In the nominal configuration (R psup =300 
cm), 1832 PMTs are placed at a radius of 3 me- 
ters in a 3-frequency icosahedral pattern, for 75% 
coverage. 



and the energy calibration parameters ao and 
ai are determined from the Monte-Carlo calcu- 
lations. 2 The energy resolution function, which 
is directly related to the total number of hits, is 
shown in Fig. [5] 

The distribution of PMT hit times (relative to 
the trigger time) can be used in the position re- 
construction algorithm to improve position reso- 
lution. For a PMT which fired at time t pm t for a 
scintillation event which occurred at time to, we 
define the time residual 

"res — Lpmt *G ™ \< ) 

where 

dpmt = vertex-PMT distance, 
n = average refractive index. 

In the absence of scattering, t res would be 
the convolution of the scintillation time spectrum 
with the PMT time resolution. Scattering will 
broaden the distribution and skew it toward later 

2 These parameters will ultimately be calibrated in the 
physical detector using radioactive sources. 
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Figure 3. Reconstructed energy (a) and position 
(b) distributions for simulated electron events. 
Position reconstruction shown using spatial infor- 
mation only. 



times. Fig. [7] shows the t res distribution defined 
by equation for simulated 100 keV electron 
events at the detector center (r = 0) and near 
the PMTs (r = R psup )- To incorporate timing in- 
formation into the reconstruction algorithm, a 2- 
dimensional probability density function (pdf) in 
t res versus event radius (shown in Fig.JSJ) has been 
derived from MC calculations, for which we can 
calculate the probability density for a given event 
with a given vertex position and time,P(x, t). A 
likelihood function is defined 

C' =R^i'P(x- f ,t r ,xl i ,t i ) (8) 

where {x},tf} are the event position and time, 



Figure 4. Position reconstruction calibration for 
simulated electron events. Position response is 
primarily dependent on the number of hit PMTs. 
For events with a small number of hits, the dark 
noise hits pull the fit toward the center. 



U} are the set of PMT positions and hit times. 
For computational efficiency, we redefine the like- 
lihood C as 



N, 



(9) 



We can now minimize equation [5] to find the 
most likely values for {x}, tf}. An x-z slice of the 
likelihood surface defined by equation is shown 
in Fig.^a) for a simulated 20 keV electron event 
at position (0, 0, R psup ), and the likelihood versus 
event time is shown in Fig. EJb), where a deep 
minimum is seen near the event position and time. 
The information in the spatial hit pattern of the 
PMTs can be added to the likelihood function 
by including to equation |^1 the position response 
using only the hit information: 



C 



NlfP(x), t r , &*o] x Ps(x- f , £) (io) 



where P s (xj,£,i) is the position response function 
shown in Fig. [3J 

Fig. 1 101 shows the reconstructed position distri- 
butions using equation 1101 and equation [3] (spa- 
tial only) for 20 keV electron events at r = R psU p- 
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Figure 5. Position resolution function derived 
from simulated electron events using spatial in- 
formation. 



The resolution for low-energy events improves by 
a factor of approximately two when the time in- 
formation is included. This is a significant im- 
provement since it implies the same background 
reduction as doubling the detector size, assuming 
20 keV background at the PMT radius. The de- 
rived position resolution functions versus kinetic 
energy for both methods are shown in Fig. ITTI 

The response functions shown in Figures El 
and ^2 will be used in the following sections to 
predict the energy and position distributions for 
signal and background events in the detector. 
The calibration requirements for energy and po- 
sition reconstruction are discussed in section [7| 

4. Solar neutrino sensitivity 

Solar neutrinos are detected by their elastic 
scattering from electrons in the target material. 
The recoiling electrons will generate scintillation 
light. The predicted differential electron recoil 
spectrum from elastic scattering of solar neutri- 
nos is 

dT e J dT , n , 

+ Kr{T')^{T' 1 T e )]dT' 
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Figure 6. Energy resolution function derived from 
simulated electron events. The energy resolution 
has little dependence on event position. 



with 

n = number of electrons in target 
\ x = e or t component of v flux 
da x 

— r = v — e cross section 
dT' 

The neutrino-electron elastic scattering cross- 
section is taken from JUj. The e and {/x, r} com- 
ponents of the solar neutrinos can be estimated 
by using the standard solar model predicted spec- 
tra and rates [20|, and the solar neutrino mixing 
parameters from a global analysis of the solar and 
reactor neutrino experiments [21] • With the pa- 
rameters from Table 0] we can calculate the pre- 
dicted solar neutrino spectrum as 

A e (T) = A(T)P ee (T) (12) 
A MT (T) - A(T)(1 - P ee (T)) (13) 

where A(T) is the undistorted solar neutrino spec- 
trum, and P ee is the electron neutrino survival 
probability at Earth calculated with the nomi- 
nal solar mixing parameters. 3 Fig. I12r a) shows 
the predicted pp solar neutrino energy spectrum. 
Fig- H2f b) shows the predicted differential elec- 
tron recoil energy spectrum for a 3 meter radius 

3 The survival probability was calculated with analysis 
code provided by the SNO collaboration. 
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Figure 7. Time residual distribution for simulated 
100 keV electron events at r=0 and the PMT ra- 
dius (r=R psup ). The difference in these timing 
distributions leads to the improved position re- 
construction. 



target of liquid neon exposed to pp neutrinos for 
1 year, and figure H^T c) shows the predicted spec- 
trum measured with the detector, for a fiducial 
volume cut of 125 cm, representing approximately 
10 fiducial tonnes of target. The total number of 
events within these cuts is 3819, which leads to 
a statistical uncertainty on the neutrino flux of 
approximately 1.6%. 



Table 4 

Solar neutrino absolute flux and mixing parame- 
ters. 

Parameter value 
pp flux 5.95 x 10 1U cm~ 2 s _1 
Am 2 7-l±S;§ x 10- 5 eV 2 

9 32 I 5±U degrees 




Figure 8. Two-dimensional time residual pdf de- 
rived from Monte-Carlo calculations used in the 
maximum likelihood position reconstruction algo- 
rithm. 
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Figure 9. (top frame) Likelihood surface in x-z 
plane for simulated 20 keV electron event at posi- 
tion (0, 0, Rpsup), and (bottom frame) likelihood 
versus event time. A deep minimum is found near 
the true event position and time. 



—140 










o"120 






Spatial vert 


=x litter 


c 
o 










ftoo 








Ttex fitter 


£ 










| 80 

To 










f£ 60 










40 




ln(a r ) = 


-0.4548 ln(T e ) + 4.7077 




20 










°C 












100 


200 300 400 


500 



T, (keV) 



Figure 11. Position resolution with spatial and 
temporal reconstruction algorithms. Including 
time information improves position resolution by 
a factor of approximately two for low energy 
events. 
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Figure 12. (a) Solar pp neutrino energy spectrum 
predicted by Bahcall et. al. (20J - (b) Predicted 
electron recoil energy distribution in liquid neon, 
(c) Predicted observed energy distribution with a 
125 cm fiducial volume cut (approx. 10 tonnes). 
Figures (b) and (c) are normalized to 1 live year. 
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5. Backgrounds 

5.1. External backgrounds 

Purification of the neon is expected to reduce 
background contamination from internal sources 
of radioactivity significantly 1 . The PMTs 
and associated hardware will contain the largest 
amount of radioactive contamination near the 
inner detector volume. To reduce these back- 
grounds we must design an active shielding region 
of liquid neon, and apply position reconstruction 
cuts to define a fiducial signal region with ac- 
ceptable background levels. With the nominal 
impurities in the PMTs, the number of decays 
from the uranium and thorium chains, and from 
40 K is large, shown in Table [S] In principle we 
would like to simulate all of the activity from 
these sources of background, and apply event re- 
construction and fiducial volume cuts to estimate 
the residual background contamination in the sig- 
nal region. This is difficult in practice given finite 
computational resources, since the total number 
of simulated events would be large. 
To decrease computational requirements, the de- 
tector response to x-ray and 7-ray photons from 
the uranium and thorium chains and from potas- 
sium decays in the PMTs is estimated by using 
the MC-calculated response functions in energy 
and radius from section 2: 

fescape(T )fldT dr 

where j(T ,r ) is the distribution of energy depo- 
sition (in kinetic energy T and radius r ) due to 
photons from the PMTs, and R(T, T' ,r,r') is the 
detector response function. f escape is a correction 
factor for photon absorption in the PMT glass, 
shown in Fig. ^3 and Vl is a solid angle correc- 
tion to account for x-ray and 7-ray coincidence 
(approximately 0.5 for x rays and 1 for 7 rays). 
We further assume 

R(T,T',r,r) = R(T,T) x R{r,r ,T') (15) 

and use the response functions shown in 
Figs. II II and IB1 4 In this scheme we need only the 

4 Note the subtle difference between R(r, r , T) and 



position and energy response functions derived, 
and a model for the 7-ray distribution at the 
PMTs due to radioactive contamination in order 
to estimate residual backgrounds. The distribu- 
tion 7(T , r ) is calculated with Monte-Carlo by 
generating the photon energy distribution from 
the uranium and thorium chains, and from potas- 
sium decays (including all x and 7 rays listed in 
the Table of radioactive isotopes [22], shown in 
Fig. I14J) and recording the resulting radial and 
energy distributions. Low-energy x-rays are most 
problematic even though they are strongly atten- 
uated by the PMT glass and the liquid neon, since 
the poor position resolution at low energies allows 
these events to misreconstruct within the fiducial 
volume. We can integrate equation ^] between 
r = and r = i-q to estimate the effective en- 
ergy spectrum due to decays in the PMTs for a 
fiducial volume cut, tq. Fig. 1151 shows these for 
several fiducial volume cuts. Fig. ED shows the 
PMT backgrounds for each contaminant, and the 
intersection with the pp energy spectrum which 
defines the analysis threshold. The achievable 
background reduction for the commercially avail- 
able PMTs would allow a low enough threshold 
(approximately 13 keV) for a feasible pp neutrino 
and dark matter experiment. 

5.2. Internal backgrounds 

The dominant internal source of background 
for the pp neutrino measurement is expected to 
be 85 Kr since it has a relatively short half- life 
(« 11 years), is a (3 emitter with energies in the 
same range (Q-value = 687 keV) as the pp neu- 
trinos, and is present in the atmosphere. Sev- 
eral other naturally-occurring radioactive con- 
taminants, expected to be less troublesome than 
85 Kr, will need to be removed from the neon in 
order to achieve acceptable background levels for 
neutrino detection, most notably uranium and 
thorium. To estimate the requirements for inter- 
nal background contamination, numbers of events 
from the uranium chain, thorium chain, and from 
85 Kr were estimated, and target concentrations 
were set which would lead to less than 1% back- 
ground to the solar pp measurement between 20- 



R(r,r',T'). 
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Figure 13. Photon escape fraction from PMT 
glass versus photon energy, for 1mm and 3mm 
thick glass. Low-energy x rays are highly attenu- 
ated by the glass. 



500 keV (the region of interest for pp neutrinos). 
These target concentrations are shown in Table 
The required concentrations for uranium and tho- 
rium, although quite low, have been achieved in 
large quantities of organic liquid scintillator [2~3*| . 
Due to the different binding energy of neon on 
carbon relative to those of most other impuri- 
ties [2U, cryogenic purification with cold traps 
is expected to be quite effective. The require- 
ments for 85 Kr are severe and ex-situ assays at 
this level will present unprecedented challenges. 
The spectral shape of j3 decays from 85 Kr, shown 
in Figure nTT a - ) can be used to extract this activity 
in-situ from the PMT data. Figure I17f shows 
the inferred statistical uncertainty on the inte- 
gral pp neutrino flux from a maximum likelihood 
(ML) separation of 85 Kr decays, 7 Be neutrino, 
and pp neutrino interactions. It is found that a 
concentration of 10~ 15 of 85 Kr can be extracted 
in-situ without significantly increasing the pp un- 
certainty, although the correlation with the 7 Be 
signal will make separation of 85 Kr from 7 Be neu- 
trino interactions difficult. 



5.3. Cosmogenics and neutron back- 
grounds 

The primary backgrounds for dark matter and 
solar neutrino experiments related to cosmic ray 
muons are fast neutrons produced in materials 
surrounding the detector and long-lived radioac- 
tivity created in the detector itself from muon 
spallation. Mei j2U has estimated the size of 
these backgrounds for a liquid neon experiment 
and concludes that they can be made negligible by 
locating the detector at sufficient depth (greater 
than 5000 meters water equivalent) to reduce the 
incoming muon flux. 

The large spherical volume of neon will also act 
as a veto for neutrons from sources external to 
the detector which penetrate a water shielding re- 
gion. Neutrons will thermalize and be captured in 
short distances in liquid neon. A fiducial volume 
cut can thus be used to remove external neutrons, 
and the coincidence signal between the scintilla- 
tion photons from neutron-neon recoils and the 
subsequent capture 7 ray can further reduce this 
background. In addition the time distribution of 
a neutron event will be distinguishable from a 
WIMP-nucleus scatter, since a neutron in liquid 
neon will scatter several times, while a WIMP will 
scatter only once. The neutron time distribution 
is thus broader than the corresponding distribu- 
tion for a WIMP-induced recoil. In this study 
the detector is assumed to be at a sufficient depth 
and have adequate shielding so that these back- 
grounds are negligible. 

6. Dark Matter Sensitivity 

In the approach taken here, the ultimate sensi- 
tivity to dark matter is determined by the limit 
on the total rate of detected dark matter inter- 
actions. This rate is a function of the WIMP- 
nucleon cross-section, the fiducial target mass, 
the analysis energy threshold, and the number of 
irreducible background events. We present here 
the predicted rates and recoil energy spectra for 
WIMP-nucleon scattering in liquid neon. 

6.1. Dark matter recoil spectra and rates 

The differential nuclear recoil spectrum for 
WIMP dark matter elastic scattering from a tar- 
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get with nuclear mass A is given by (in events 
per unit energy per unit time per unit detector 
mass) [2SJ 



with 



dR 
dQ 



with 



Po 

Vo 

m 



x 



niN 

HQ) 

T{Q) 



;F 2 (Q)T(Q) 



(16) 



= unquenched nuclear recoil 

kinetic energy 
= WIMP-nucleus ES cross section 
= local dark matter density 
= speed of Sun around galactic center 
= WIMP mass 

= reduced WIMP-nucleus mass 
m x mN 
m x + rriN 
— Nuclear mass 
= form factor 



erf- 



-Vq X 



erf- 



v v 
minimum WIMP velocity 
earth's velocity. 



We use here the standard nuclear form factor 



HQ) 



exp --(gs) 

qr n 2 



(17) 



with the numerical values given by |27) : 



m p = proton mass 



(19) 



where we consider only the spin-independent in- 
teraction. 

Figure 1181 shows the nuclear recoil spectrum 
calculated with eauationll6lfor the input param- 
eters defined in Table [7| for several target nuclei. 
A potential advantage over a xenon-based detec- 
tor is the lower sensitivity to the achievable en- 
ergy threshold due to the shape of the recoil spec- 
trum. Figure El shows the ratio of expected dark 
matter rates above threshold for xenon and neon 
targets, where it is seen that the relative xenon 
sensitivity drops quickly as the energy threshold 
is increased. The nuclear form factor correction 
(Eqn. I17|l tends to suppress the predicted event 
rate for large A nuclei with increasing kinetic en- 
ergy, so that for experimentally achievable en- 
ergy thresholds, much of the gain obtained by ex- 
ploiting the coherence predicted by equation 1181 
for large- A target nuclei is lost. With a 10 keV 
threshold, for example, the additional sensitiv- 
ity with xenon is only approximately a factor of 
6, which can be overcome with neon by using a 
larger target mass. 



Table 7 

Parameters used 



for calculation WIMP recoil 



spectrum and rates. 


Parameter 


Value 


Po 


0.3 GeV c~' z 


vo 


220 km s" 1 



qr n = 6.92 x 10- 3 A 1/2 Q 1/2 (1.UA 1/3 ) 
qs = 6.92 x 10" 3 A 1/2 Q 1/2 0.9, 

jl is spherical Bessel function, and Q is the nu- 
clear recoil energy in keV. 

To compare experiments with different target 
nuclei, we express the nuclear cross-section in 
terms of the nucleon cross-section [2H1 



(T = Cp 



1 + ?71 X /TOJV 

1 + m x I m p 



A 2 



(18) 



6.2. Distinguishing nuclear and electron 
recoils in neon 

Assuming we have designed the experiment so 
that the backgrounds to the pp neutrinos after 
threshold and fiducial volume cuts are small, the 
dominant background to potential dark matter 
events will arise from the pp neutrinos themselves. 
In this case the ultimate sensitivity to dark mat- 
ter will be limited by the expected dark matter 
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interaction rate and the separability of dark mat- 
ter events from neutrino events. The PMT time 
distribution can in principle be used to separate 
these classes of events, since the intrinsic scintil- 
lation time distribution is different for electron 
and nuclear recoils (see Table |2Jl - The PMT time 
distribution predicted by the Monte-Carlo for 20 
keV electron events and 80 keV nuclear recoil 
events from WIMP elastic scattering is shown in 
Fig. I20f a). These events would yield similar ef- 
fective energy distributions due to the quenching 
of the scintillation light from nuclear recoils. To 
separate these classes of events, the ratio f prompt 
of prompt to total hits is calculated. These dis- 
tributions are shown in Fig. I20f b~) . The effect of 
PMT noise is to decrease the prompt fraction. 

We can now use these distributions to statisti- 
cally separate nuclear recoils from electron events. 
The separation efficiency was estimated by gen- 
erating 1000 MC datasets, each with 3000 solar 
neutrino events, corresponding to approximately 
one live year. We define the extended negative 
likelihood function 

N cvcntB 

£ = -2(N e + N x ) J2 WePe(f prompt )+ (2Q) 
^X^xii 'prompt)] 

with P e and P x the probability distributions for 
f prompt shown in Fig. 120T b) . We can then mini- 
mize equation [3UJ to find the most likely values of 
iV e and N x . Since equation[2UJis an extended like- 
lihood, the s — a parameter uncertainties can be 
calculated by finding the iso-likelihood contours 
which satisfy 

and include the sampling uncertainties in our 
dataset, so that they include the uncertainties on 
deduced incident particle fluxes. 

Equation [2U] was minimized (with the ROOT 
TMinuit package [3UJ) for each of the 1000 simu- 
lated datasets. The fitted numbers of electron 
and WIMP events were found to peak at the 
input values, providing a check of signal analy- 
sis. The 90% CL upper limit on the number of 
WIMP events in the data set is estimated from 



this technique to be approximately 1.3 events, 
which provides background discrimination of ap- 
proximately 0.9965. We can use this upper limit 
along with the predicted nuclear recoil rate from 
WIMP scattering (which is the convolution of 
equation 1161 with the energy response function 
of Fig. [S] corrected for 25% quenching) to plot 
the sensitivity contour in WIMP-nucleus cross- 
section versus WIMP mass. This is shown in 
Fig. ESI along with current and projected ex- 
perimental results. The limit corresponds to 
a WIMP-nucleon spin-independent cross section 
(<Tp) of -10- 46 cm 2 for a 100 GeV WIMP mass. 
Details on the selection of optimal analysis win- 
dows and cuts are provided in the following sec- 
tion. 

6.3. Optimizing detector sensitivity 

The optimal energy threshold and fiducial vol- 
ume for the dark matter and pp experiment will 
in general be different, since the energy spectra of 
these two event classes differ. To simultaneously 
optimize the experiment for dark matter and pp 
neutrinos, an optimum fiducial volume and en- 
ergy threshold are found for each by maximiz- 
ing the total predicted event rate above thresh- 
old as a function of fiducial volume cut. It was 
found that an energy window of 15 keV above 
threshold provided the best WIMP sensitivity for 
a 100 GeV WIMP, and this window was employed 
throughout the analysis. Fig.l21lshows the deter- 
mination of the optimal fiducial volume for neu- 
trinos and dark matter for a 300 cm detector. It is 
found that the sensitivity for dark matter is opti- 
mized by selecting a smaller fiducial volume and 
thus a lower analysis threshold than for the so- 
lar neutrino analysis, since the energy spectrum 
of the WIMP interactions is much harder. The 
optimized thresholds and fiducial volumes, which 
are different for the dark matter and pp compo- 
nents of the experiment, have been used in the 
analysis. From Fig. I21f b). the relatively large 
region over which the dark matter sensitivity is 
optimized is due to a trade-off between a lower 
analysis threshold (for smaller fiducial radii) and 
better WIMP-neutrino separation with increasing 
analysis threshold (for larger fiducial radii). This 
effect can be used a systematics check on the data 
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from the liquid neon experiment and allows some 
flexibility in ultimately determining the analysis 
cuts. The optimized sensitivities for the pp and 
dark matter experiments are shown versus detec- 
tor size in Fig. [221 From this figure a 300 cm ra- 
dius detector is found to provide good sensitivity 
for both measurements. Little additional sensi- 
tivity is found for detector sizes much larger than 
400 cm in radii, since the measurements become 
limited by position resolution. 
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Figure 14. 7 and x-rays from uranium (a) and 
thorium (b) chain decays, and from potassium 
(c). Shown are the number of events per keV per 
chain decay [in (a) and (b)], and per decay [in (c)]. 
Approximately 5% of each chain decay results in 
a photon in the range 0-50 keV. 
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Figure 15. External background from contamina- 
tion with D737KB PMTs at 300 cm, for different 
fiducial volume cuts. 
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Figure 17. In-situ determination of 85 Kr back- 
ground, (top) T c ff spectra predicted for pp and 
7 Be neutrinos, and for 85 Kr. (bottom) Uncer- 
tainty on pp flux from in-situ determination ver- 
sus Kr concentration. 



Figure 16. Backgrounds from uranium, thorium, 
and potassium for a PMT sphere with radius 300 
cm, and a 125 cm fiducial volume. 
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Figure 18. Nuclear recoil spectra (unquenched ki- 
netic energy) from WIMP elastic scattering on Xe 
(A=131), Ge (A=73) and Ne (A=20) calculated 
with equation^] Rates have been calculated as- 
suming dp = 1CP 42 cm~ 2 and m x = 100 GeV, 
and are in units of keV _1 kg _1 day _1 . The points 
are from a calculation by the CDMS collabora- 
tion [2HI 
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Figure 19. (a) Relative rates above threshold per 
unit mass of material for a 100 GeV WIMP on Xe 
and Ne. Contours show relative rates above the 
effective energy threshold, assuming 25% quench- 
ing for Xe and Ne. (b) Relative rate above 
threshold assuming the same quenching (25%) 
and threshold for both Xe and Ne. The relative 
rate decreases rapidly with increasing threshold, 
and is approximately 6 for a 10 keV threshold. 
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Figure 20. Distinguishing electron versus nu- 
clear recoil events based on PMT time distribu- 
tion, (top) PMT time distribution for 10000 sim- 
ulated electron recoil and nuclear recoil events, 
(bottom) Distribution of f prompt for these sets of 
events. 



Figure 21. Optimal fiducial volume determina- 
tion for the pp neutrino (a) and dark matter (b) 
experiments for a nominal PMT sphere radius of 
300 cm, and 1 year of livetime. 
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Figure 23. Limits on the WIMP-nucleon spin- 
independent cross-section versus WIMP mass for 
current and planned experiments and projected 
sensitivity for liquid neon. Data from current ex- 
periments and projected XENON sensitivity from 
Gaitskell and Mandic pTT] , 



Figure 22. Sensitivity for (a) the pp neutrino 
and (b) dark matter experiments versus radius 
of PMT sphere (R psU p). The fiducial volume and 
thresholds are optimized for each value of R psU p- 
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Table 5 

Decay rates in D737KB PMTs. 



Chain 


Specific activity 


Total mass in PMTs 


Decays/year 




(Bq/g) 


(g) 




U 


178.9 x 10 3 


0.0553 


3.122 x 10 11 


Th 


39.14 x 10 3 


0.0553 


6.83 x 10 10 


K 


30.30 


110.52 


1.057 x 10 11 



Table 6 

Internal radioactivity requirements. The target mass and concentrations are for a background of less 
than 1% of the pp solar neutrino rate with an energy threshold of 20 keV. 



Nuclide 


Activity 


Isotopic 


Specific 


Target 


Target 




source 


abundance 


activity 


mass 


Cone. 








(Bq/g) 


( n §) 


(g/g) 


Kr 


8b Kr 


1.5e-ll 


1.55 x 10 b 


3.7 x 10" 4 


1 xl0~ ib 


Th 


232 Th chain 


1.00 


39.14 x 10 3 


3.1 


2.2 xl0~ 17 


U 


235 U, 238 U chains 


0.993, 0.007 


178.9 x 10 3 


8.5 


6.2 xl0~ 17 
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7. Detector calibration requirements 

Assuming radioactive backgrounds to the pp 
flux measurement can be reduced such that they 
contribute less than 1% uncertainty, the domi- 
nant systematic uncertainties will be due to un- 
certainty of the acceptance of the solar neutrino 
signal. These uncertainties will arise through un- 
certainty of the energy threshold and fiducial tar- 
get mass. Table |H1 outlines a nominal error bud- 
get assuming a target pp flux uncertainty of 1%, 
where the dominant systematic uncertainties are 
assumed to each contribute an equal 0.35% to the 
total uncertainty. The dominant fiducial mass 
(volume) uncertainty arises from uncertainty in 
the difference between the true and reconstructed 
radii of events. The fractional radial uncertainty 
€r — =j^- is related to the fiducial target uncer- 
tainty 4^ by 



e« = 



,/av; 



1 - 1. 



(22) 



The requirement of a 2.5 mm position uncertainty 
for a 200 cm fiducial volume will require a cali- 
bration source which can be positioned to that ac- 
curacy, an order of magnitude improvement over 
that currently achieved by SNO and similar ex- 
periments. This requirement could be relaxed if 
the overall uncertainty on the pp flux is allowed 
to increase, if other systematic uncertainties are 
lower than shown in Table [5J or if the statisti- 
cal uncertainty is reduced by increasing detector 
livetime. 

The uncertainty due to the absolute energy 
scale calibration can be evaluated by calculating 
the change in area of the spectrum of Fig. I12f b) 
when the threshold changes by one standard er- 
ror. The minimum threshold uncertainty versus 
threshold is shown in Fig. [21 for the target pp 
flux systematic uncertainty due to energy scale 
calibration of 0.35%. Shown in Table |S] are the 
requirements for calibration of the absolute en- 
ergy scale (near threshold) and the position re- 
construction offset for a resultant 0.35% uncer- 
tainty on the pp flux, assuming a 35 keV thresh- 
old. It is interesting to note that a large scintil- 
lation yield versus neon temperature dependence 
[approx (8%/K)] is found in 0, which leads to the 



requirement that the neon be temperature stable 
at the level of AT< 125 mK to satistify the re- 
quirements of Tabic [9] 

The calibration requirements for the dark mat- 
ter component of the experiment are much less 
stringent. Fiducial volume uncertainties of a few 
percent achieved by current generation experi- 
ments are adequate. An uncertainty of 20-30% 
in the energy scale calibration for nuclear recoil 
events near the threshold (~ 10 keV) leads to a 
factor of two change in the deduced WIMP cross- 
section, and will be achievable with calibration by 
a neutron source. 
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Figure 24. Absolute energy scale calibration un- 
certainty (as a fraction of the threshold energy) 
versus threshold for a 0.35% systematic uncer- 
tainty on the measured pp flux. 



8. Design criteria 

Table EH summarizes measurements which 
must be made to allow for an optimized design 
of a liquid neon pp neutrino and dark matter ex- 
periment. An "X" in either the dark matter or 
pp neutrino column denotes the measurement pri- 
marily affects this component of the experiment. 
Dependence of the ultimate detector sensitivity 
on these parameters is discussed in the following 
sections. 
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Table 10 

Requirements for design criteria of a liquid neon detector for both solar neutrinos and dark matter. 
Response characteristics marked with "X" need to be measured (or shown to be within range) before the 
detector design can be optimized for either the dark matter or solar neutrino experiment. 



Response 


Dark matter Solar v 


e~ scint. yield 


X 


X 


e~ scint. prompt/late ratio 


X 




Scintillation yield linearity 


X 


X 


nucleus scint. yield 


X 




nucleus scint. prompt /late ratio 


X 




late WS time constant 


X 




Rayleigh scattering 






absorption length 


X 


X 


neon temperature stability 




X 


PMT long-term stability 


X 


X 


PMT QE 


X 


X 


PMT threshold 


X 


X 


PMT dark noise 


X 


X 



Table 8 

Error budget for pp flux measurement with 



R D 



=300 cm, 1 year livetime. 



Source 


Relative uncert. 




(%) 


statistics 


0.7 


energy scale 


0.35 


vertex accuracy 


0.35 


uranium background 


0.35 


thorium background 


0.35 


krypton background 


0.35 


total 


1.0 



Table 9 

Absolute energy scale and position calibration re- 
quirements for 0.35% contributions to the pp flux 
systematic uncertainty, assuming a threshold of 
3 5 keV. 

energy scale (AE/E) 

vertex accuracy (AR/R) 

vertex accuracy (at R=200 cm) 



< 1% 

< 0.12% 
2.5 mm 



8.1. Photon absorption in liquid neon 

Absorption of the primary scintillation photons 
in the liquid neon will both reduce the total light 
output and introduce systematic effects related 
to event position and direction. To estimate the 
effects of absorption on total light output, a set 
of Monte-Carlo calculations were performed for a 
range of photon absorption lengths. Fig.l25lshows 
the light output (relative to no absorption) for 
scintillation events at the center of a 300 cm ra- 
dius detector for absorption lengths ranging from 
100 to 10 5 cm. The effect of the 60 cm Rayleigh 



scattering length is to greatly increase the optical 
path length for photons, and thus decrease the 
tolerable amount of absorption. To obtain a rel- 
ative light yield of 95%, an absorption length of 
approximately 300 m is required. 

The absorption coefficient (a) for scintillation 
photons can be estimated by adding the absorp- 
tion coefficients for the various absorbers present 
in the liquid: 



23 




Absorption length (cm) 



Figure 25. Relative light output (normalized to 
no absorption) versus absorption length for events 
at the detector center. 



These levels are moderate compared to that re- 
quired for 85 Kr, and most impurities are expected 
to be easily removed from the neon through the 
use of cold traps. The problem of hydrogen dif- 
fusion into the neon is expected to be mitigated 
by the very low diffusion rate at liquid neon tem- 
perature. 

Table 11 

Total absorption cross sections for 80 nm photons 
and impurity limits for a 300 m absorption length. 
Absorber Cross section at 80 nm ji 



(xlQ- 18 cm 2 ) (xlQ- 12 ) 

H 2 10.5 2.2 

H 2 21 9.8 

N 2 12 26.8 

2 17 21.6 



where 



Mi 

P 

N A 

H 



Molar mass of absorber 
liquid neon density 
Avogadro's constant 
Mass fraction of absorber 
absorption cross section 



If we assume the absorption is equally divided 
among N absorbing materials 5 , we can set the 
mass fraction limit (7^) for each absorber: 



M, 



NpN A ai\ 



(24) 



Neon has no significant absorption for photons 
with wavelengths above the first ionization limit 
at 575 nm, and so is transparent to its scintil- 
lation light A comprehensive list of photo 
absorption cross-sections has been compiled by 
Gallagher et. al. [33]. Shown in Table ITU is a list 
of dominant absorbers present in neon, along with 
the cross-sections for photo absorption of 80 nm 
photons from 3*3*j , and the required impurity lim- 
its for a total photon absorption length of 300 m. 

5 This is likely an oversimplification since some materials 
will pose more of a problem than others, but this nonethe- 
less provides an estimate of the tolerable contaminations. 



8.2. Dependence on optical and detector 
parameters 

Many of the parameters shown in Ta- 
bles H an d El are either estimates inferred from 
other measurements or have not been quanti- 
fied with systematic uncertainties. The sensi- 
tivity of the solar pp neutrino measurement to 
these parameters is quite modest. To study 
the sensitivity of the dark matter experiment to 
these inputs, a set of Monte-Carlo calculations 
were performed by varying the parameters which 
will affect the detector performance by a nomi- 
nal amount (±25% in most cases, or a physical 
bound). Figure [5^1 shows the deduced sensitivity 
for a 100 GeV WIMP versus detector size, for the 
ranges of optical and detector parameters of inter- 
est. This provides a benchmark for the required 
accuracy with which these parameters need to 
be measured, and points out the parameters to 
which we are most sensitive. The prompt to total 
light yield for electrons and nuclear recoil events 
provide separation between these two classes of 
events, and are found to be the dominant factors 
in determining the ultimate WIMP sensitivity. It 
is interesting to note that the WIMP sensitivity 
varies only weakly with total light yield, nuclear 
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recoil quenching factor, and the Rayleigh scatter- 
ing. This is somewhat fortuitous since measuring 
these absolute quantities experimentally in the 
laboratory with high precision is difficult. Ex- 
periments are currently underway at Los Alamos 
National Laboratory and Yale University to mea- 
sure these properties. 

Figure |23 shows the effect of changing the total 
PMT radioactivity by ± one order of magnitude 
from the nominal values for PMTs with commer- 
cially available ultra-low background glass. This 
is found to also have a relatively small effect on 
the estimated WIMP sensitivity. 
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Figure 26. Dependence of WIMP sensitivity on optical and detector parameters. Central values are 
from model simulations and are shown as dashed curves. Shown are the ranges in sensitivity for (a) 
electron prompt to total light ratio varied by ±25%, (b) nuclear recoil prompt to total light ratio varied 
by +4.8,-25%, (c) PMT noise rates of and 1000 Hz, (d) Light quenching factor for nuclear recoils varied 
by ±25%, (e) total scintillation yield varied by ±25%, and (f) Rayleigh scattering varied by ±25%. The 
WIMP sensitivity depends most critically on the prompt to total light ratio for electrons and nuclear 
recoils. 
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Figure 27. WIMP sensitivity for an order of mag 
nitude change in the PMT backgrounds. 



9. Discussion 

In summary, we have demonstrated the pos- 
sibility for a simultaneous dark matter and pp 
neutrino experiment using liquid neon, assum- 
ing nominal scintillation characteristics and back- 
ground contamination levels with detailed Monte- 
Carlo calculations. This approach to direct dark 
matter detection has the advantage that WIMP 
event discrimination is based only on PMT sig- 
nals, allowing for a simple detector design which 
can in principle be scaled to a very large target 
mass. This detector design also has the potential 
advantage that external sources of background 
are now well-separated from the inner target ma- 
terial, and can be dealt with through position 
reconstruction. Assuming light attenuation and 
internal source radioactive backgrounds can be 
controlled, the large target mass achievable with 
neon may lead to the best sensitivity of any of the 
direct dark matter searches. We have shown that 
form factor suppression for high-A target nuclei 
leads to less of an advantage over lower-A target 
nuclei than that expected by the cross-section co- 
herence alone. A specific list of issues which must 
be addressed in order to determine the ultimate 
detector sensitivity has been provided. An ex- 



perimental program is underway to address these 
issues. 
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